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Abstract —This paper deals with the problem of adaptive 
multidimensional/multichannel signal detection in homogeneous 
Gaussian disturbance with unknown covariance matrix and 
structured (unknown) deterministic interference. The aforemen¬ 
tioned problem extends the well-known Generalized Multivariate 
Analysis of Variance (GMANOVA) tackled in the open literature. 
In a companion paper, we have obtained the Maximal Invariant 
Statistic (MIS) for the problem under consideration, as an 
enabling tool for the design of suitable detectors which possess the 
Constant False-Alarm Rate (CFAR) property. Herein, we focus on 
the development of several theoretically-founded detectors for the 
problem under consideration. First, all the considered detectors 
are shown to be function of the MIS, thus proving their CFARness 
property. Secondly, coincidence or statistical equivalence among 
some of them in such a general signal model is proved. Thirdly, 
strong connections to well-known simpler scenarios found in 
adaptive detection literature are established. Finally, simulation 
results are provided for a comparison of the proposed receivers. 

Index Terms —Adaptive Radar Detection, CFAR, Invari¬ 
ance Theory, Maximal Invariants, Double-subspace Model, 
GMANOVA, Coherent Interference. 


I. Introduction 
A. Motivation and Related Works 

T HE PROBLEM of adaptive detection has been object of 
great interest in the last decades. Many works appeared 
in the open literature, dealing with the design and performance 
analysis of suitable detectors in several specific settings (see 
for instance 12 and references therein). 

As introduced in a companion paper, herein we focus on 
a signal model which generalizes that of GMANOVA 0 by 
considering an additional unknown double-subspace structured 
deterministic interference. Such model is here denoted as I- 
GMANOVA. The I-GMANOVA model is very general and 
comprises many adaptive detection setups as special instances, 
ranging from point-like targets (resp. interference) 0 to 
extended ones 0, from a single-steering assumption to a 
vector subspace one 0, (6), and the GMANOVA model itself 
0, only to mention a few examples. We recall that attractive 
modifications of GMANOVA have also appeared in the recent 
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literature 0, 0, focusing on the design of computationally- 
efficient approximate ML estimators when the unknown signal 
matrix is constrained to be diagonal 0 or block-diagonal 0. 

In the case of composite hypothesis testing, the three 
widely-used design criteria are the Generalized Likelihood 
Ratio Test (GLRT), the Rao test, and the Wald test 0. Their 
use is well-established in the context of adaptive detection 
literature 0, ED-IJ21; more important they are known to 
share the same asymptotic performance 0. However, in the 
finite-sample case their performance differ and their relative 
assessment depends on the specific hypothesis testing model 
being considered. Such statement holds true unless some 
specific instances occur, such as in fl4l . where it is proved that 
they are statistically equivalent in the case of point-like targets 
and a partially-homogeneous scenario. Other than the afore¬ 
mentioned detectors, in the context of radar adaptive detection 
it is also customary to consider their two-step variations, 
with the two-step GLRT (2S-GLRT) being the most common. 
Those are typically obtained by designing the detector under 
the assumption of the a known disturbance covariance matrix 
and replacing it with a sample estimate based on the so-called 
secondary (or signal-free) data OH. 

Furthermore, a few interesting alternative detectors for com¬ 
posite hypothesis testing are the so-called Durbin (naive) test 
lfl6l and the Terrell (Gradient) test ITT) . These detectors have 
been shown to be asymptotically efficient as the aforemen¬ 
tioned well-known criteria. The same rationale applies to the 
Lawley-Hotelling (LH) test liT8Tl . Though these detectors are 
well-known in the statistics field, the development and appli¬ 
cation of these decision rules is less frequently encountered in 
radar adaptive detection literature , e.g., 021, ED. The reason 
is that an important prerequisite for a wide-spread application 
of an adaptive detection algorithm consists in showing its 
CFARness with respect to the nuisance parameters; in this 
respect, the assessment of such property in radar adaptive 
detection literature has been somewhat lacking. 

Of course, the use of GMANOVA model in the context 
adaptive radar detection is not new and dates back to the the 
milestone study in (13), where the development and analysis 
of the GLRT was first proposed. A similar work was then 
presented years later in (20), where the focus was on the 
design of a compression matrix prior to the detection process, 
aimed at reducing the computational burden and minimizing 
the performance loss with respect to the standard processing. 
More recently, GLRT, Rao and Wald tests were developed and 
compared under the GMANOVA model 122 . along with some 
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other heuristic detectors. Unfortunately, albeit the CFARness 
of the proposed detectors was proved, no clear connection to 
the MIS was established. More importantly, no (structured) 
deterministic interference was considered in all the aforemen¬ 
tioned works; the proposed I-GMANOVA model is aimed at 
filling such gap. 

We point out that the closest study to ours (in terms of 
interference accounting) is the work in CD, where range- 
spread and vector subspace targets and interference are con¬ 
sidered; however the sole GLRT and 2S-GLRT are derived 
and analyzed. Similarly, a Rao test (and its two-step version) 
is recently obtained in Il22l . It is worth noticing that the model 
considered by the aforementioned works is included in the I- 
GMANOVA model, and can be readily obtained by assuming 
a canonical form for the right-subspace matrix of both the 
signal and the interference. 

Summarizing, in our opinion the huge (but scattered) lit¬ 
erature on adaptive detection in many case-specific signal 
models and the presence of several detectors (developed on 
theoretically-solid assumptions) for generic composite hypoth¬ 
esis testing problems lacks a comprehensive and systematic 
analysis. First, such analysis may help the generic designer in 
readily obtaining a plethora of suitable adaptive detectors in 
some relevant scenarios which can be fitted into the considered 
I-GMANOVA model. Secondly, the development of detectors 
closed-form expressions under I-GMANOVA model may al¬ 
low to easily claim some general statistical equivalence results 
than those already noticed in some special instances (see e.g., 
El, El, ED). Thirdly, the available explicit expression for 
each detector allows for a systematic analysis of its (possible) 
CFARness (under a quite general signal model). The latter 
study is greatly simplified when knowledge of the explicit 
form of the MIS is available for the considered problem; in 
this respect, the derivation of the MIS and its analysis, object 
of a companion paper, fulfills this need. 


B. Summary of the contributions and Paper Organization 

The main contributions of the second part of this work 
are thus related to detectors development and CFAR property 
analysis and can summarized as follows: 

• Starting from the canonical form obtained in our com¬ 
panion paper, for the general model under investigation 
we derive closed-form expressions for the (i) GLRT, (ii) 
Rao test, (Hi) Wald test, ( iv ) Gradient test ( v ) Durbin 
test, ( vi ) two-step GLRT (2S-GLRT), and (mi) LH test. 
As an interesting byproduct of our derivation, we show 
that Durbin test is statistically equivalent to the Rao 
test for the considered (adaptive) detection problem, thus 
extending the findings in fl9l . obtained for the simpler 
case of a point-like target without interference. Similarly, 
we demonstrate the statistical equivalence between Wald 
test and 2S-GLRT, thus extending the works in El 
and Il23l . concerning the special instances of point-like 
targets (no interference) and multidimensional signals, 
respectively; 

• The general expressions of the receivers are exploited 
to analyze special cases of interest, such as: (a) vector 


subspace detection of point-like targets (with possible 
structured interference) a, 0, ei, (b) multidimen¬ 
sional signals GDI, ED, (c) range-spread (viz. extended) 
targets J4), l24l . (25], and (d) standard GMANOVA (i.e., 
without structured interference) ED, El- In such special 
instances, possible coincidence or statistical equivalence 
is investigated among the considered detectors; 

• Exploiting the matrix pair form of the MIS obtained in 
part one, we show that all the considered detector can 
be expressed as a function of the MIS, thus proving their 
CFARness with respect to both the covariance of the dis¬ 
turbance and the deterministic (structured) interference; 

• Finally, a simulation results section is provided to com¬ 
pare the proposed detectors in terms of the relevant 
parameters and underline the common trends among 
them. 

The remainder of the paper is organized as follows: in Sec. EI 
we describe the hypothesis testing problem under investiga¬ 
tion; in Sec. m we obtain the general expressions for the 
detectors considered in this paper and we express them as a 
function of the MIS; in Sec. [IV] we particularize the obtained 
expressions to the aforementioned special instances of adap¬ 
tive detection problems; finally, in Sec. [V] we compare the 
obtained detectors through simulation results and in Sec. m 
we draw some concluding remarks and indicate future research 
directions. Proofs and derivations are confined to an additional 
document containing supplemental materiaQ 


II. Problem Formulation 

In a companion paper, we have shown that the considered 
problem admits an equivalent (but simpler) formulation by 
exploiting the so-called “canonical form”, that is: 



Z — A |_S £ 0 

Z = ABeC 


0(Vf xr 

-N 


C + N 


( 1 ) 


where we have assumed that a data matrix Z E C NxK 
has been collected. Also, we have adopted the following 
definitions: 


1 Notation - Lower-case (resp. Upper-case) bold letters denote vectors (resp. 
matrices), with a n (resp. A n ,m ) representing the n-th (resp. the (n, m)-th) 
element of the vector a (resp. matrix A); R N , C N , and W Nx ^ are the 
sets of TV-dimensional vectors of real numbers, of complex numbers, and 
of TV x TV Hermitian matrices, respectively; upper-case calligraphic letters 
and braces denote finite sets; E{-}, (-) T , (-)t, Tr [•], ||-||, 9^{-} and S{}, 
denote expectation, transpose, Hermitian, matrix trace. Euclidean norm, real 
part, and imaginary part operators, respectively; 0 nxM (resp. In) denotes 
the NxM null (resp. identity) matrix; On (resp. In) denotes the null (resp. 
ones) column vector of length TV; vec(TVf) stacks the first to the last column 
of the matrix T \I one under another to form a long vector; det(A) and 11 A| | f 
denote the determinant and Frobenius norm of matrix A; A<g)B indicates the 
Kronecker product between A and B matrices; denotes the gradient of 

scalar valued function fix) w.r.t. vector x arranged in a column vector, while 
i ts transpose (i.e. a row vector); the symbol means “distributed as”; 
x ~ CA/*iv(/x, X) denotes a complex (proper) Gaussian-distributed vector x 
with mean vector \a E C Nxl and covariance matrix X E C NxN ; X ~ 
CNnxm{A, B , C) denotes a complex (proper) Gaussian-distributed matrix 
X with mean A E <C NxM and Cov[vec(A)] = B <S) C\ Pa denotes 
the orthogonal projection of the full-column-rank matrix A, that is Pa = 
[A(A^ A) -1 At], while P^ its complement, that is P ^ = (I — Pa)- 
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. A = [E t E r ] G C NxJ , where E t 4 [l t 0 tx(iv _ t) ] T 

and E r = [0 rxt L are the (known) left 

subspaces of the interference and useful signal, respec¬ 
tively (we have denoted J = r + t); 

• B s = [Bj 1 B t ] T , where B t .i G C txM and B G 
C r x A1 are the (unknown) interference (under Hf) and 
useful signal matrices, respectively; 

• C = [Im o Mx(K-M)\ G C MxK is the (known) right 
subspace matrix associated to both signal and interference 
in canonical form; 

• TV is a disturbance matrix such that N ~ 

CJ\f nxk(Onxk,Ik,R), where R G C NxN is an 

(unknown) positive definite covariance matrix G5). 

We recall that the detection problem in (|T]i is tantamount to 
testing the null hypothesis B = 0 rX M (viz. ||B|| f = 0, 
denoted with Hq) against the alternative that B is unrestricted 
(viz. ||.B|| f > 0, denoted with Hi), along with the set of 
nuisance parameters B f l and R. 

In the present manuscript we will consider decision rules 
which declare Hi (resp. Ha) if <h(Z) > V (resp. <1>(Z) < 

rj), where <!>(•) G C NxK —> R indicates the generic form 

of a statistic processing the received data Z and r) denotes 
the threshold to be set in order to achieve a predetermined 
probability of false alarm ( Pf a ). 

As a preliminary step towards the derivation of suitable 
detectors for the problem at hand, we also give the following 
auxiliary definitions: 

• bn G R rMxl and bj G R rMxl are obtained as bji = 
Ji{6} and bj = A{ 6 }, respectively, where we have 
defined b = vec (B)-, 

• 0 r = [ b^, bj ] G R 2rMxl is the vector collecting 
the parameters of interest; 

• 9 S = [O^a Oj.b\ T S ]^(2tM+Af 2 )xt j s t jj e yec _ 

tor of nuisance parameters containing: (a) 9 sa = 
[ bjj j bfj ] T G K 2tMxl where bt : R and bt.i are the 
vectors obtained as b t _R = R { b t } and b t j = 3{6 f }, 
respectively, where b t = vec {BA (i.e. B t i under Hi)\ 
(b) contains in a given ordei[j the real and imaginary 

parts of the off-diagonal entries together with the diagonal 
elements of i?; 

. 0 = [ Oj 9j ] T G R( 2JM + Ar2 ) xl is the overall 


unknown parameter vector; 

^ r iT ^ 

• 9o = 9j (j 9j 0 , with 9 S fi denoting the Maximum 

Likelihood (ML) estimate of 9 S under Hq and 9 r o = 

0 2 r M (that is, the true value of 9, under Hq)', 

^ r ,t ^ ^ 

• 9i = 9j- { 9j { , with 9, \ and 9 s l denoting the 

ML estimates of 9 r and 9 S , respectively, under Hi. 

The probability density function (pdf) of Z, when the hypoth¬ 
esis Hi is in force, is denoted with /i(-) and it is given in 
closed form as: 


fi{Z-B s ,R) =TT- NK det(R)~ K 
x exp (—Tr [R~ 1 (Z — AB S C)(Z — AB S C)^]) , (2) 


while the corresponding pdf under Hq, denoted in the follow¬ 
ing with /o(-)> is similarly obtained when replacing AB S C 
with E t BtfiC in Eq. 0. In the following, in order for 
our analysis to apply, we will assume that the condition 
(K — M) > N holds. Such condition is typically satisfied 
in practical adaptive detection setups US. 


A. MIS for the considered problem 

In what follows, we recall the MIS for the hypothesis testing 
under investigation, obtained in our companion paper. The 
mentioned statistic will be exploited in Sec. [Til] to ascertain 
the CFARness of each considered detector. Before proceeding 
further, let 


Y,i = 


Im 

0(K-M)xM 


v c ,± 


0 Mx(K-M) 

Ik-m 


(3) 


and observe that P c t = (V^iV^) and P^ t = (V^V^)- 
Given these definitions, we denote (as in Part I): (i) Z, = 
(ZT C , 1 ) G C NxM , (H) Z c , x 4 (ZV C , 2 ) G C JVx ^- M ) and 
(fit) S c ± (Z Ct± zl ± ) = (ZP^Zt) G C NxN . 

It has been shown in our companion paper that the MIS is 
given by: 


T(Z C ,S C ) 


rp A f 7 ! C-l 

J-a ~ |^ 2.3 ^ 2.3 

T b 4 {z] ^ 

z 2 f Sn 1 Z 2 



}1 


J < N 


(4) 


J = N 


where Z 2 . 3 = (Z 2 - S 23 S^ Z :i ) and S 2 . 3 = {S 22 - 
SnS^Sn). Also, we have exploited the following parti¬ 
tioning for matrices Z c and S c : 


'Zl' 


's 11 

S 12 

S l 3 

Z 2 

, S c = 

S 21 

S 22 

s 23 

Z 3 


S:n 

S 32 

S 33 


where Z x G C txM , Z 2 G C rxM , and Z 3 G C^ N ~ J '> xM , 
respectively. Furthermore, S V) , (i, j) G {1, 2, 3} x {1, 2, 3}, is 
a sub-matrix whose dimensions can be obtained replacing 1, 
2 and 3 with t, r and (N — J), respectively Additionally, for 
notational convenience, we also give the following definitions 
that will be used throughout the manuscript: 


^23 = 


Z 2 

Z 3 



( 6 ) 


Finally we recall that, for the I-GMANOVA model, the in¬ 
duced maximal invariant equals T p = B^ R 2 ! B G C MxM , 
where J? 2 3 is analogously defined as S 2 . 3 when S r is replaced 
with the true covariance R. 


III. Detectors design 

In this section we will consider several decision statistics 
designed according to well-founded design criteria. Initially, 
we will concentrate on the derivation of the well-known GLRT 
(including its two-step version), Rao and Wald tests ©. 


2 More specifically, 0 s b = where S(-) denotes the one-to-one 

mapping providing 0 S ^ from R. 


3 Hereinafter, in the case J = N, the “3-components” are no longer present 
in the partitioning. 
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Then, we will devise the explicit form of recently used 
detection statistics, such as the Gradient (Terrell) test EH, 
the Durbin (naive) test which have been shown to be 
asymptotically distributed as the three aforementioned detec¬ 
tors (under very mild conditions). Finally, for the sake of 
completeness, we will obtain the LH test for the problem at 
hand, following the lead of fl5l . 


sKN 


-K 


A. GLR 

The generic form of the GLR in terms of the complex¬ 
valued unknowns is given by |9l : 

ma *{b s m} fi(Z ; B S ,R ) 
max{ Bl ^ H } fo(Z: B t ,o,R) 

First, it can be readily shown that the ML estimate of R under 
Hi (resp. under Ho), parametrized by B s (resp. B t ,o) is: 

Ri(B s ) = K-'iZ- AB S C)(Z- AB S C/ ( 8 ) 
Ro(B t , 0 ) = K~ l (Z - Ek B t>0 C)(Z - E t B tfi C)* (9) 

After substitution of Eqs. ([ 8 ]) and <[9]) in /i(-) and fo(-), 
respectively, the concentrated likelihoods are expressed as: 

f 1 (Z-B s ,R l {B s )) = {K/{ Tte)) 
xdet [(Z- AB S C)(Z~ AB S C)/ “ (10) 

fo(Z-,B t>0 ,Ro{B ti0 )) = (K/( ne)) KN 
xdet[(Z-E t B tfi C)(Z-E t B tfi C)l] ~ K (11) 

Then the ML estimates of B s and B t o under Hi and Ho, 
respectively, are obtained as the solutions to the following 
optimization problems: 

B s 4 argmindet[(Z - AB S C)(Z - AB S C )/; (12) 

B s 

B t o 4 argmindet[(Z - E t B tfi C)(Z - E t B tfi C )t]. 

Bt,0 

(13) 

It has been shown in jl5l that the optimizers have the closed 
form: 

B s = (A t S'- 1 A)- 1 A t S'- 1 ZC t (CC t )- 1 ; (14) 

Bt ,o = (E\ S' 1 Et)-'E\ S~ l ZC\C&)-\ (15) 

Substituting Eqs. (IT4l) and < 1 15 b into ( ITOb and i | lib , respectively, 
provides (after lengthy manipulations): 

fi(Z; B a ,Ri) = (K/(ne)) KN det [S C ]~ K 

x det [I M + (Z wl V cA ) f P^ (Z W1 V cA )) ~ K (16) 

/o(Z; B ti0 ,Ro) = (K/(ne)) KN det[S c ]-* 

x det [I M + (Z wl V Ct i) f Pi 0 (Z W1 V cA )j~ K (17) 

where we have defined Ai = (S c 1 ^ 2 A), Ao = (S c 1 ^ 2 E t ) 
and Z W1 = (S c * 1 Z), respectively. Finally, substituting 
Eqs. dT 6 ]) and (fl7T> into Eq. ([TJ 1 and after taking the fc-th root, 
the following explicit statistic is obtained: 

det [Im + (Z wl V cA ) t PUZ wl V c ,i)\ 


tglr — 


det [Im + ( Z W 1 V C>1 )^ P ^( Z wl V cA )\ 
det [Ik + Z^ W1 P^ o Z wl P c t] 
det[J^ + Z \ V1 Pf Z W 1 P C t] 


( 18 ) 

( 19 ) 


where the last expression follows from Sylvester’s determinant 
theorem |26l . Furthermore, we observe that Eq. ( IT8l ) can 
be also re-arranged in the following useful equivalent forms 
(again obtained via Sylvester’s determinant theorem): 

tgir = det [Im - D~ 1/2 (Z wl V cA /V A (Z wl V cA )Df ll2 }~ 1 

( 20 ) 

= det[J m + D- 1/2 (Z wl V c ^V A (Z wl V c ,i)D^ 2 ) (21) 

where P A = (Pa 1 - Pa 0 ) and A = [I M + 
(Z W1 V C , i)l P 4 . (Z Wi Vc,i)\, respectively. Finally, it is worth 
noticing that Eq. (ITH i is in the well-known Wilks’ Lambda 
statistic form E2). Moreover, the latter expression generalizes 
the GLR in fL5ll to the interference scenario (i.e., t ^ 0). 

For the sake of completeness, we also report the closed form 
ML estimates of R obtained under Ho and Hi (after back- 
substitution of (fril l and (IT5l > in Eqs. © and ©, respectively): 

Ri(B s ) = K- l [S c + (Z- S}/ 2 P Al Z W1 )P C] 

x(Z~S 1 J 2 P Al Z wl )/, (22) 

Ro(B t ,o) = K~ l [S c + (Z- S l J 2 P Ao Z W1 )P & 

x(Z ~S l J 2 P Ao Z wl )/, (23) 

and underline that we will use the short-hand notation R, in 
what follows. Finally, before proceeding further, we state some 
useful properties of ML covariance estimates (later exploited 
in this paper) in the form of the following lemma. 

Lemma 1. The ML estimates of R under Hi and Ho satisfy 
the following equalities: 

Rf 1 A = K Sf 2 A = K [Sf 1 E t Sf [ E r ] (24) 

R a 1 E t =KSf 1 E t (25) 


Proof: Provided as supplementary material. ■ 

CFARness of GLRT: Using the expression in Eq. (1 1 8b . we 
here verify that t g \ T can be expressed in terms of the MIS (cf. 
Eq. (I}), thus proving its CFARness. Indeed, it can be shown 
thaf 


( Z wl V cA / Pi 0 ( Z wl V c y ) = z \ 3 Sf \ Z 2 . 3 
( Z wl V c > 1 ) 1 p £ 1 ( Z wl V c , i ) = zl S3-3 1 z 3 , 


+ z \ Sf / Z 3 , 
(26) 

(27) 


from which it follows 


det [Im + T a + Tf\ 
det [Im + T&] 


(28) 


which demonstrates invariance of the GLR(T) with respect to 
the nuisance parameters and thus ensures CFAR property. 


B. Rao statistic 

The generic form for the Rao statistic is given by S): 


dlnfi(Z-d) 


del 


[I- 1 (do)] 


e r ,e r 


a in fi ( Z - d ) 


0=00 


80 r 


(29) 


0=00 


4 The proof of the aforementioned equalities is non-trivial and thus provided 
as supplementary material. 
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where 


I{0) = E 


j^n/ 1 (Z;0)dln/ 1 (Z;0) 


(30) 


86 dO T 

denotes the Fisher Information Matrix (FIM) and 
[I~ 1 (6)\e r ,g r indicates the sub-matrix obtained by selecting 
from the FIM inverse only the elements corresponding to the 
vector 6 r . It is shown (the proof is provided as supplementary 
material) that the aforementioned statistic is given in closed 
form as: 

n- 1 E r f° e! R- 1 Z d 


tr 


= Tr 


Z d 0 Ro Er ^22 E r P0 Z d ,0 P(7t 


where we have partitioned r° = (A^ R 0 1 A) 1 


r° = 


11 


12 


(31) 


(32) 


and, r°-, (i,j) S {1,2} x {1,2}, is a sub-matrix whose 
dimensions can be obtained replacing 1 and 2 with t and r, 
respectively. Additionally, we have defined: 

z d ,0 =(z- sy 2 P Ao S - 1 ' 2 Z P c t ^ 

Eq. ( OTb can be rewritten in a more familiar (and convenient) 
way. Indeed, it is provec@ that: 

Pet, (34) 


(33) 


R- 1/2 Z dA P c ,=Pi o Z 


where Aq = (i? 0 
and also 


d ,0 Pct 
E t ) and Z wo = (R Q 


Z), respectively. 

Pa o *0 1/2 E r r 22 El Po 1/2 Pi = ( Pa , -PAo), (35) 

where A\ = R 0 1 ^ 2 A and Aq = R (j 1 ^ 2 E t , respectively. 
Therefore, an alternative form of / rao is obtained substituting 
Eq. ( PH into Eq. (OH and exploiting (PH . thus leading to the 
compact expression: 

(36) 


tr. 


— Tr \Z\ V0 (P Al P Ao ) Z wo P c 


Ct 


CFARness of Rao Test: We now express the Rao statistic as 
a function of the MIS, aiming at showing its CFARness. To 
this end, we first notice that Eq. (OH can be rewritten as: 


t — Tr 

‘'rao — -L 1 

{z wo v c y\Pi^ Pi)(z wo v c y 


(37) 

Moreover, exploiting the equal itief] 




,.i Piyz wo v c ,f) = k {43 1 

^23 


-ZlSf 1 

z 23 {i m + z \ 3 si z 23 )~ l zy si 

Z 23 

f. 




(38) 

(Z W 0 V cA y Pi (Z W 0 V Ct i) = K [z\ Si 

Z 3 


-z\ si 

zyi M + zl si za^zI si z 3 

f. 

(39) 


Eq. 


can be rewritten as: 

frao = K Tr [T a — ( T a + T b ) 


x (Jm + T a + T b ) 1 ( T a + Tb ) 

+ T b (I M +T b )- 1 T b } (40) 

which is only function of the MIS, thus proving its CFARness. 


5 The proof is provided as supplementary material. 

6 Their proof is provided as supplementary material for this manuscript. 


C. Wald statistic 

The generic form for the Wald statistic is given by (9): 

(0r,i - e rfi ) T {[I~ l {6{)\e r ,er}- 1 (»r,l - 0r,o) • (41) 


It is shown (the proof is provided as supplementary material) 
that the aforementioned statistic is given in closed form as: 


fwald = Tr 


'z^KPiS-^Er V\ 2 E}S^ 2 P^ 0 Z w1 


Pet 

(42) 


where FT indicates the (i,})-th sub-matrix of T 1 = 
(A.1 Rf 1 A)~ l , obtained using the same partitioning as in 
Eq. (OH for matrix r°. The above expression can be rewritten 
in a more compact way, as shown in what follows. Indeed, the 
inner matrix in Eq. (PH is rewritten a^|: 

K Pi 0 S c “ 1/2 Er f 22 E}Sf l / 2 Pi 0 = V A (43) 

which then gives: 

iwaid = Tr [Z^ V1 P a Z W1 Pet] (44) 


CFARness of Wald Test: Finally we prove CFARness of 
Wald statistic. First, it is apparent that Eq. (PH can be rewritten 
as: 


Uaid = Tr [(Z wl V cA y{Pi o - Pij Z W1 V C>1 ] . (45) 

Secondly, exploiting Eqs. (OH and (OH (as in the GLR case), 
Eq. (PH is rewritten as: 


fwald — Tr 


Sfl Z 2 .3 


J 2.3 ^2.3 


= Tr [Ta 


(46) 


Therefore f wa id depends on the data matrix uniquely through 
the MIS (actually, only through the first component). 


D. Gradient statistic 

The Gradient (Terrell) test requires the evaluation of the 
following statistic llT7l . (28): 


dhi h{Z;0) 


dOT 


(0r, 1 0 r , o) 


e=e o 


(47) 


The appeal of Eq. (PH arises from the fact that it does not 
require neither to invert the FIM nor to evaluate a compressed 
likelihood function under both hypotheses (as opposed to 
GLR, Wald, and Rao statistics). As a consequence, this formal 
simplicity can make the Gradient statistic easy to compute. 
Moreover, under some mild technical conditions, such test is 
asymptotically equivalent to the GLR, Rao and Wald statistics 

Gll- 

It is shown (the proof is provided as supplementary material) 
that the Gradient statistic is given in closed form as: 


^grad 


4 5R |Tr [zyKPyS^^Er f\ 2 E}.Rf' Z d , o Pet] } 


( 48 ) 


7 The proof is provided as supplementary material. 
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where Z d , o is given in Eq. (l33l >. The expression in Eq. ( l48l > 
can be cast in a more compact form, as shown below. First, 
we notice that the following equality holdfj 

Ro 1 z d , o Pc t = S-^P^ ( Si ' 2 Rq 1 ) z P c t (49) 

which, after substitution in Eq. (l48l ) and exploitation of 
Eq. (01, gives the final form 

fgrad = K{Tr [Z' W1 V A (Si/ 2 R- 1/2 ) Z wo P c t]}, (50) 


where identical steps as for Wald test have been exploited. 

CFARness of Gradient Test: First, Eq. (l50l) can be readily 
rewritten as: 


fgrad = 5R {Tr [(Z wl V Ctl ) t (Pi 0 - Pjj 
R^ 1/2 ) Z W 0 V C ,^} . 

Moreover, exploiting the following equal itie{3 


((Z wl V cA y Pi 0 S l J 2 R/ 1/2 Z W0 V c a) 

= K (T a + T b )[I M - (Im +T a + + T b )] 

((Z wl V cA y P^ Sl' 2 Pq : l/2 ZwoV cA ) 

= KT b [I M - (Im + T a + Tf,) _1 (T a + T 6 )] 

it thus follows that: 


(51) 


(52) 

(53) 


igrad = K {Tr [JlT a (I M - (Im + T a + Tb)" 1 

x(T a + Tb))]} (54) 

which shows that also the Gradient test satisfies the CFAR 
property. 


E. Durbin statistic 

The Durbin test (also referred to as “Naive test”) consists 
in the evaluation of the following decision statistic ESI: 



where the estimate 0 r ,oi is defined as: 

0r,oi = argmax/i(Z;0 r ,0 SiO ). (56) 

O r 

In general, the Durbin statistic is asymptotically equivalent to 
GLR, Rao and Wald statistics, as shown in ESI- However, for 
the considered problem, a stronger result holds with respect 
to the Rao statistic, as stated by the following theorem. 

Theorem 2. The Durbin statistic for the hypothesis testing 
model considered in Eq. 0 is statistically equivalent to the 
Rao statistic. Therefore, the test is also CFAR. 

Proof: Provided as supplementary material. ■ 

It is worth noticing that the present result generalizes 
the statistical equivalence observed between Rao and Durbin 
statistics for the simpler scenario of point-like targets and 
single-steering assumption in lfl9l . On the other hand, Thm.[2] 
proves that such result holds for the (very general) hypothesis 
testing problem considered in this work. 

s The proof is deferred to supplementary material. 
lJ The proof is provided as supplementary material. 


F. Two-step GLR (2S-GLR) 

It is also worth considering a two-step GLR (2S-GLR), 
which first consists in evaluating the GLR statistic under 
the assumption that P is known and then plugging-in a 
reasonable estimate of P. The GLR statistic for known P 
can be expressed in implicit form as ES): 

max Bs /i (Z;B S ,R) 
max Bt 0 f 0 (Z;B t: o,R) ' 

The ML estimates of B s and B t 0 are more easily obtained 
from optimizing the logarithm of /i(-) and / 0 (-), respectively, 
that is: 


-K ln(n N det [P]) 

— Tr[P~ 1 (Z’ — AB S C)(Z — AB S C)/ (58) 

— K ln(7r w det[P]) 

- Tr[P _1 (Z - E t P t>0 C)(Z - E t P t>0 C) t] (59) 

Maximization of Eqs. (l58l > and ( l59l > with respect to B s and 
P tj o, respectively, can be obtained following the same steps 
employed in ITT31 and thus it is omitted for brevity. Therefore, 
after optimization, the following statistic is obtained (as the 
logarithm of Eq. (l57ll): 

Tr[Z+ P" 1/2 (P Ai - P Ao ) R-^ZPct] (60) 

where we have defined A\ = (P _1 ^ 2 A) and Ao = 

(P -1 / 2 P t ), respectively. We recall that the expression in 
Eq. (l60l) depends on P. We now turn our attention on 
finding an estimate for the covariance P. Clearly, in order 
to obtain a meaningful estimate to be be plugged in both the 
numerator and denominator of Eq. (l57l >. such estimate should 
be based only on signal-free data (also commonly denoted as 
“secondary data”). 

It is not difficult to show that the covariance estimate based 
only on secondary data is given b)0: 

R sd = (K-M)~ l S c . (61) 

Thus, substitution of Eq. (|6TT) into Eq. (l60l > leads to the final 
form of 2S-GLR: 


Tr[Z t \JK - MSf 1 ' 2 Va Sf 1/2 s/K - M Z P c 


c tj 


^2s— glr — Tr 


Zf S ~ 1/2 Va Sf 1/2 ZP ( 


Ct 


oc 

(62) 


From direct comparison of Eqs. 6H» and (l62l >. a general 
equivalence result is obtained, stated in the form of the 
following lemma. 


Lemma 3. The 2S-GLR statistic is statistically equivalent to 
the Wald statistic. Therefore, the test is also CFAR. 


The aforementioned lemma extends the statistical equiva¬ 
lence observed between 2S-GLR and Wald statistics in the 
simpler cases of point-like targets ED, range-spread targets 
fl24l and multidimensional signals ll23l . 


10 It should be noted that the same result would be obtained by considering 
(i.e., the ML estimate under l~i\) as the signal-free covariance estimate. 
Indeed, since Eq. 4601 depends only on R through the quantities (R~ 1 A) 
and (R~ 1 Et) an thus Lem. |T| could be exploited to obtain the same final 
statistic. 
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G. Lawley-Hotelling (LH) statistic 

Finally, for the sake of a complete comparison, we also 
consider (and generalize) the simpler statistic proposed in 031 
pag. 37] as a reasonable approximation to GLR. Indeed, Wilks’ 
Lambda form of GLR in Eq. (|2T1 > can be rewritten as: 

tgk = det[7 m + D; 1/2 {BliA'S-'A) B 3 

- Bl^Els-'Ei) B t , 0 } £>r 1/2 ] (63) 

where we exploited (CC^)' 1 = (CC^) -1 / 2 and closed- 
form estimates for B s and B'\ 0 (cf. Eqs. ( 1 1 4] and (fl5l >. 
respectively). As the number of samples I\ grows large, we 
can invoke approximation S c « (K — M)R, that is, the 
sample covariance based on secondary data will accurately 
approximate the true covariance matrix. Accordingly, we can 
safely approximate 


D i ~ Im , 

B s ps (A^R- 1 A)- 1 A t J?- 1 C t (CC t )- 1 , 

B t ,o ps {ElR^Et)- 1 E t R' 1 C t (CC t )" 1 , (64) 

since (A'I'S'- 1 A) ps (K - M)~ 1 (A^R~ 1 A) and 
(EjS^Et) ss (K — M)~ 1 {E\R~ 1 E t ), respectively. 
Therefore, based on these approximations, it is apparent that 
the second contribution within the determinant in Eq. (|63| > will 
be a vanishing term as the number of observations increases. 
Hence, GLR statistic will be given by the determinant of an 
identity matrix plus a small perturbing term. 

Additionally, we remark that when T £ H MxM j s a 
small perturbing matrix, det [Im + Y] can be (accurately) 
approximated at first order as nf=ii(l ps 1 + v i = 
1 + Tr[Y], where v t denotes the z-th eigenvalue of Y. Based 
on those reasons, we formulate the LH statistic as: 


tih — Tr 


d; 1/2 (z wl v cA y v A (z wl v Ctl ) d~ x/2 

= Tr [(Z wl V cA )i V A (Z wl V cA ) Df 1 ] 


(65) 


CFARness of LH statistic: The CFARness is proved by 
using Eqs. (l26l > and (l27l > from Sec. IIII-Al within Eq. (l65l ). thus 
obtaining: 

t lh = Tr[T a {I M + T b )~ 1 ] . (66) 

Finally, in Tab. |T| it is shown a recap table, summarizing all the 
considered detectors and their respective expressions in terms 
of the MIS in Eq. @. 


IV. Detectors in Special Cases 
A. Adaptive (Vector Subspace) Detection of a Point-like Target 

In the present case we start from general formulation in 
Eq. (U} and assume that: (i) t = 0 (i.e., there is no interference, 
thus J = r and A = \l r 0 rx (7v-r)] T G <C Nxr ); ( ii ) M = 1, 
i.e., the matrix B collapses to a vector b £ C Jxl and (in) 
c = [l 0 • • • 0] £ C lxK (i.e., a row vector). Such case 

has been extensively dealt in adaptive detection literature 0, 
0 , 113, |J29), ED. The hypothesis testing in canonical form 
is then: 


Clearly, since in this case M = 1 holds, (K — 1) vector 
components are assumed signal-free, that is, Z admits the 
partitioning Z = [z p Z s ] = [z c Z c ,_l] , where z p de¬ 
notes the signal vector related to the cell under test and the 
columns of Z s represent the secondary (training) data. Also, 
Pa 0 = 0 nxn (resp. P^ Q = In) holds, because of the 
absence of the structured interference. In the latter case, it 
can be shown that the simplified projector form holds: 


Pet 


1 

Or-i 


0 


T 

K -1 




( 68 ) 


Given the results in Eq. (l68l >. it can be shown that S c = 
Z S Z\ and R a = -^S Q , where S a = (.z p z | + Z S Z\) hold, 
respectively. In some cases we will also use the Sherman- 
Woodbury formula ||26l applied to Sq 1 , that is: 


S 


-l 

o 


= S7 1 - 


S c 1 z p z^S 
1 A Zp S c 



(69) 


GLR: In the specific case of M = 1, the following form of 
the GLR is obtained from Eq. (l20l) : 


fglr — 


1-77 


A zlSf 1/2 P Al Sl 

V = — - 


- 1 / 2 . 


1 + Zp Sc 


(70) 


since we have exploited Dq -a do = (1 + z p iZ Pi i) and 
(Z wl V cA V Va (Z wl V c p) -A (.zl A P Al z p p ), where z p>1 = 
(Sf 1 / 2 z p ). Clearly, t g i r is an increasing function of rj, the 
latter thus being an equivalent form of the statistic and 
coinciding with the so-called multi-rank signal model GLR 
described in 0, EE). 

Rao/Durbin statistic: For the present scenario, Eq. ( |36| ) 
specializes into: 


t rao = Tr[Zl 0 P Ai Z W 0 P c ,} = Tr(zl(R 0 1/2 P Ai Rf 1 / 2 )z p } 
cx z+ Sf 1 A(A^ 0 - 1 A)- 1 AtS' 0 -iz p 4 Vraa (71) 

Eq. (1711) can be further simplified by exploiting the Woodbury 
identity in ( |69) (and similar steps as in 11131 1). thus obtaining 
the following simplified form of the Rao statistic: 


? 7 rao — 


1 

i ■+■ z p,i z p,i 


Z p,lPAi z p,l 
1 “f Z p pP Al Zp i 


(72) 


Finally, for r = 1 (i.e., a single-steering case) (A —> a £ 
C Arxl ), Eq. (l72l ) reduces to: 


Iz/S 1- 

1 « 2 /( atS c lfl ) 

[1 + zlS^Zp] 

1 1 -tc-n i4sc _i ar 

1 + Z p Oc Z p ( atS -i a ) 


(73) 


which coincides with the well-known Rao statistic for the 
single-steering case developed in E3. 

Wald/2S-GLR statistic: Starting from Eq. ( I44| i. we particu¬ 
larize the Wald statistic as follows: 


fwaid = r Pc{Zl vl P Al Z W iP c i i ] = Ti[z p l P Al z Pt i] (74) 
= zlS^AiA^S^A^A^S^Zp 


H 0 : Z = N 
Hi'. Z = Abc + N 


( 67 ) 
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Table I 

Detectors comparison and their functional dependence of the MIS (viz. CFARness). Auxiliary definitions : T a+b = (T a + T b ) and 

Di 4 [i M + (Z W1 V C>1 )t P z . (Z W1 V C} i)]. 


Detector 

Standard Expression 

MIS function 

GLR 

det [.Do]/ det[£)i] 

det [Im + T’a+b]/ det [Im + T b \ 

Rao/Durbin 

Tr[Z^o [L’a 1 ~ ^A 0 ) Ao Pet ] 

K Tr [T a — T a+b (I M + Ta+b) 1 Ta+b + Tb{lM + Tb)~ 1 Tb] 

Wald/ 2 S-GLR 

r ^ r [^Wl(^ > A 1 — Pa 0 ) Z W1 P c t] 

Tr [T a \ 

Gradient 

»{Tr [Z f wl (P Al - P Ao ) (S c 1/2 Rf 1/2 ) Z wo P ct ]} 

Sf {Tr [KT a {Im - ( I M + T' a+ 6 )- 1 T a+6 )] } 

LH 

Tr[(Z wl V Cll )t (P Al - Pa 0 ) (Z wl V cA ) D- 1 ] 

Tr [T a (I M + T,,)- 1 ] . 


In the special case r = 1 (i.e., a single-steering case), Eq. (El 
becomes: 


fwald — 


4 s c 1q I 

a^Sc 1 a 


igrad = 3 ?|Tr 

= ^{ z l 1 P Al (S 1 J 2 R- 1/2 )z p , 0 } 


= A'5R(4 5- 1 A(A t S'- 1 A)- 1 A t 5' 0 ” 1 2 p } 


where z p $ = 
exploiting Eq. 


— 1 / 2 

(#0 7 %)• It is interesting to note that. 


tgrad 


= if K 


= if 


the gradient statistic is rewritten as: 

f (AtS-tA )- 1 (Atg-tz p ) 

\ 1 + z\ Sc'zp 

(A+S-tzp) 


l + ziSc b 


(77) 


is given as: 


(75) 

n 0 : 

Z = A 

bio 

°r_ 

1 

Ri : 

Z = A 

blp 

b T 


c+AT 
c+lV 


(78) 


Filter (AMF) ifTTI. Oil. 

Gradient statistic: In this case the gradient statistic in 
Eq. ( l50l > specializes into: 

Zl.PA^S^R^ZwoPct]} (76) 


Clearly, since in this case M = 1 holds, (if — 1) vector 
components are assumed signal-free, that is, Z admits the 
partitioning Z = \z p Z s ~\ = [z c Z Ci ±J, where z p de¬ 
notes the signal vector related to the cell of interest and the 
columns of Z s represent the secondary (or training) data. 
In the latter case, it can be shown that the same simplified 
projector form in Eq. (l68l i holds. Given the results in Eq. (1681) 
, it can be shown that S c = Z s Z\ and Rq = j^So, 
where S 0 = (S c +Sl /2 Pji 0 Sf 1/2 z p z+S , “ 1/2 P^sl' 2 ) hold, 
respectively. In some cases we will also use the Sherman- 


Woodbury formula B261 applied to S 0 and consider the 


product S 0 z p , which provides: 


3-1 


-’0 


Zrt - 


S-'z v -{S-^Pi o S-^)z p 

<D -L < 

Ao k 


(79) 


ziS~ 1/2 PLSc 1/2 z v 


1 I z t a~ 1 / 2 p± e-!/2 
l-l -ZpOc Z P 


where in last line we have omitted 5i{-} since Eq. ( 1771 is 
formed by Hermitian quadratic forms (at both numerator and 
denominator); thus it is always real-valued. Therefore Eq. (1771 
is statistically equivalent to Kelly’s GLR in Eq. (ITOl . 

LH statistic: We recall that for point-like targets, the 
condition M = 1 holds. Therefore the LH test is statistically 
equivalent to the GLRT since the operators Tr[-] and det[-] are 
non-influential when applied to a scalar value. This follows 
since the expressions in Eqs. (|2T1 > and ( l65l are thus related by 
a monotone transformation. 


B. Adaptive Vector Subspace Detection with Structured Inter¬ 
ference 

In the present case we start from general formulation in 
Eq. ([U and assume that: (i) M = 1, i.e., the matrices B 
and B t i collapse to the vectors b € C rxl and b t i £ C txl , 
respectively; (ii) c = [l 0 0] £ (C 1 ^ (i.e., a 

row vector). Such case has been dealt in ED- Given the 
aforementioned assumptions, the problem in canonical form 


GLR: In the specific case of M = 1, the following form of 
the GLRT is obtained from Eq. (l20l : 


Air — 


1 


1 - 7 ?’ 


A Z p,l (-PAi Pa 0 ) Z p p 

V = -5—"i-> (80) 

I + z p,i Pa 0 z pA 


since we have exploited D 0 —> do = (1 + 

zlSf 1/2 Pi o Sf 1/2 z p ) and (Z W1 V c ,i)t P A (Z W1 V cA ) ->• 
z l.\{ p M -Pa 0 ) z p, 1 ’ where z P ,i - (Sf 1/2 z p ). Clearly, Eq. 
( l80t is an increasing function of rj, which can be thus seen as 
an equivalent form of the GLR. 

Rao/Durbin statistic: For the present scenario, Eq. 
specializes into: 


Aao — Z p o (Pax P Ao ) z P ,C 


( 81 ) 


where z pfi = ( Rf 1,2 z p ). 

Wald/2S-GLR statistic: Starting from Eq. (l44l >. we particu¬ 
larize the Wald statistic as follows: 


fwaid — Z p,i (Pa l Paq) z 


■p, 1 


(82) 
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Gradient statistic: In this case the gradient statistic in 
Eq. (150} specializes into: 

tgrad = 5R {4, ( P Al - P Ao ) (Sl /2 R- [/2 ) z p ,o} (83) 
We now rewrite Eq. (183} as: 

tgrad = KK{zl S~ 1 / 2 (P Al - P Ao ) (5 c 1/2 So ') z p } (84) 

Exploiting Eq. (|79] > and observing that (P Al ~ Pa 0 )Pa 0 = 
Pa i — Pa 0 holds, Eq. (184} is expressed as: 


^grad — P 5ft ' 


zlS, 


- 1/2 


(P Ai -Pao)S, 


- 1/2 


i + 4 s~ lf 2 p|„ s'” 172 


= a: 


4 s c 1/2 (p Al - p A o ) s~ c 


1/2 


i + 4 s c -1/2 


pi- 0-1/2 


(85) 


where in last line we have omitted 5ft{ -} since Eq. (185} is 
formed by Hermitian quadratic forms (at both numerator and 
denominator); thus it is always real-valued. Therefore Eq. ([85} 
is statistically equivalent to GLR in Eq. (180b . 

LH statistic: As in the case of no-interference in Sec. II V-Al 
the condition M = 1 holds. Therefore the LH statistic is 
statistically equivalent to the GLR. 


C. Multidimensional Signals 


In the present case we start from formulation in Eq. ([[} and 
assume that: (*) t = 0 (i.e. there is no interference, meaning 
J = r), (ii) A = E r = I N (thus J = r = N) and (in) 
C = [Im 0 mx(k-m)\- Such case has been dealt in ifTOl . 
l23l . Thus, the hypothesis testing in canonical form is given 
by: 


Ho : Z = N 
H t : Z = BC + N 


( 86 ) 


Clearly, since in this case J = N holds, (K — M ) vector 
components are assumed signal-free, that is, Z admits the 
partitioning Z = [. Zm Z s J = [Z c Z 0) jJ, where Zm 
denotes the signal matrix collecting the cells containing the 
useful signals and the columns of Z s are the training data. In 
the latter case, it can be shown that the simplified projector 
form holds: 


Pet = 




0 


Mx(K-M) 


0(K-M)xM ®(K-M)x(K-M)_ 


(87) 


Given the results in Eq. (168} . it can be shown that S, — Z s Z\ 
and Rq = 4So, where Sq = (ZmZ^ m +Z s Z\) holds, respec¬ 
tively. Also, it is not difficult to show that P A[ = P Al = In 
and P A = P A = Onxn, respectively. 

GLR: In order to specialize GLR expression we start from 
Eq. ( IT8} . Indeed, it can be easily shown that: 


det[/M + {Z W iV c a)^ {Z W iV c a)\ 

~ det [I M + (Z wl V cA )lP^(Z w1 V c a)} 
= det[/M + l Z M ] 

= det [I M + S-^Zm Z^ 5" 1/2 ] 

= det [S c + ZmZIj] / det [S'c] 


where we have exploited Pj^ = Onxn and Sylvester’s 
determinant theorem in third and fourth lines, respectively. It 
is apparent that the latter expressions coincide with those in 
IfTOl Eqs. (18) and (20)], respectively. 

Rao/Durbin statistic: For the present setup Eq. (136} spe¬ 
cializes into: 

frao = r I I [Zl V0 PA 1 Z W0 P c -f] = Tl' [Z^R q 1 Z Pet] 

= KTr[(ZP c ,ySo 1 (ZP c ,))=KTr[Zl I So 1 ZM} (89) 

which coincides with the specific result obtained in l23l , which 
was originally derived as a modified two-step GLRT procedure 
in IfTOl . 

Wald/2S-GLR statistic: Starting from Eq. (l44l >. we particu¬ 
larize the Wald statistic as follows: 


fwaid — — Tr [Z^S C 1 Z Pet] 

= TrKZPc^S-^ZPct)] = TtIzIs^Zm] (90) 


which coincides with the specific result obtained in 11231 . 

Gradient statistic: In this case the gradient statistic in 
Eq. (Boll specializes into: 


tgrad = 5ft {Tr [z; i P Al (5 c 1/2 Po 1/2 )Z wo P Ct ] } 

= K 5ft {Tr [Z^Sq 1 Z P ct ] } = K Tr 




(91) 


It is interesting to observe that in this specific scenario, 
Gradient statistic coincides with Rao statistic in Eq. (1891) . 

LH statistic: In this specific instance, LH statistic in 
Eq. (165} specializes into: 


<M 


t\h — Tr (Z wx V^i) x 

(Zw 1 V c ,i)'P£ 1 {Z wl V c , 1 j)~ JI ] = Tr f Z f M S~ 


(92) 

Zm 


since P A = 0 ; y x ..y (resp. P A| = In) for multidimensional 
signal setup. From inspection of the last line, it is apparent that 
LH statistic coincides with Wald/2S-GLR statistic in Eq. (190} 
for this specific scenario. 


D. Range-spread Targets 

In the present case we start from general formulation in 
Eq. 0 and assume that: (i)t = 0 (i.e., there is no interference, 
thus J = r); (ii) r = 1, thus the matrices A and B collapse 
to a = [1 0 ••• 0] T G C Nxl and b £ C lxM (i.e., a 

row vector), respectively; (Hi) C = [Im 0 mxK-m\- Such 
case has been dealt in j4j, l24l . Il25l . Therefore, the hypothesis 
testing in canonical form is given by: 


Po : Z = N 

Pi : Z = abC+N ' 


(93) 


Additionally, (K — M) vector components are assumed signal- 
free, that is, Z admits the partitioning Z = [ Z, Z s ] where 
Z, £ C v x M comprises the cells containing the extended 
target and Z s £ (£ N x( K - M ) collects the secondary (training) 
data. In the latter case, the following simplified projector form 
holds: 


Pet = 


Im 0 Mx(K-M) 

0(K-M)xM 0 (K-M)xM 


( 88 ) 


(94) 
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Based on the structure of Eq. ( l94l ). it follows that S, = Z S Z\ 
and R 0 = j^So, where So = ( Z e Z e + Z S Z\). Moreover, it 
can be shown that P ai and P a( (where we have analogously 
defined a± = (S c 1 ^ 2 a) and di = (R 0 1 ^ 2 a)) assumes the 
following simplified expression: 

„ S~ 1 / 2 aa^S - 1/2 „ S- 1 / 2 aa ^~ 1/2 

a^S^a '■ ' 

(95) 

In some cases, we will use the Woodbury identity applied to 
S~\ that is: 

So 1 = s - 1 - s~ l z e ( I M + z\s~ x z e y l z\ s~ x . m 


GLR: Aiming at particularizing the expression of the GLR 
for the present case, we follow the same derivation as in fl5l to 
show that Eq. ( l20l > can be specialized exploiting the equalities 


{Z wl V cA y P ai (Z wl V cA ) = 


(ztsya)(z$sy a y 


(a^S c 1 a) 

Do = Im + Z\S~ l Z e , 


(97) 

(98) 


thus obtaining f g i r = [ 1/(1 — 77 )], where: 

a {^S~ x Z e ) [I M + Z\syZ e ]^ ( Z\S~ x a) 

V ~ (at S^a) ' ^ 

The result in Eq. (199b is obtained after substitution of Eqs. d23 
and (l98b into Eq. (l20b and exploiting Sylvester’s determinant 
theorem. Such GLR fomFI (as t g \ T is a monotone function of 
77 ) corresponds to that found in |[T5l . 

Rao/Durbin statistic: In the present case Eq. (f36b special¬ 
izes into: 


Oao — P r [Zl V0 P ai Z W0 P(jt} 

= Tr[(ZP c ^R~ 1/2 P iil R- 1/2 (ZP c ,)\ 

= Tr[Z e t R-^P^R-' x ' 2 Z e \ 

KTx[Z\ sy a ysyz e \ = K \\Z\syaf_ 

(at5' ( / 1 a) (atS' ( / 1 a) 


( 100 ) 


Eq. (11 OOl i is recognized as the result found in 124) . 

Wald/2S-GLR statistic: Starting from Eq. (l44b . we particu¬ 
larize the test as follows: 


fwald — 1t:[Z\ vl P ai Z wl P c \] 

= Tr[(ZP c ^Sy/ 2 P ai Sy/ 2 (ZP c i)] 

= Tr[Z e t S-V 2 P^S-V* z e ] 

_ T,[ZtS- l aysyZ e \ _ WZtsyaf 

(a^Sc 1 ^ (aiS'(T 1 a) 

which agrees with the result in |j24 | and can be shown to 
coincide with the generalized AMF proposed in a thus 
extending the theoretical findings in [fill . 


11 It is worth pointing out that an alternative (equivalent) form of GLR 
was obtained in a, gd for the range-spread case. The aforementioned 
expression can be simply obtained starting from general formula in Eq. 
CD, straightforward application of Sylvester’s determinant theorem and 
exploitation of the simplified assumptions of range-spread scenario. 


Gradient statistic: In this case Eq. (l50l > reduces to: 


fgrad 


= 3t {Tr [zt, 1 P ai (S c 1 / 2 fi 0 ” 1/2 )Z wo P c t] } 


= K 


= K 


3? {Tr [Z\ S- 1 aa^Sy l Z e ] } 
(atS' t T 1 a) 

Kt{[zisya]* [Ztsya]} 


( 102 ) 


(atS/ 1 a) 

where we have used 5ft{Tr[mnt]} = SR [min], with m and 
n being two column vectors of proper size. Furthermore, by 
exploiting Eq. (l96b . the following equality holds 


[z\ S^a) = [I M + Z\ SyZej-^Zi S~ x a) (103) 
which, substituted into Eq. (1 1 02b . gives: 

[ZtS- x a]'\I M + Z\SyZ e \- x [Z\S- l a] 


^grad — 


(atS'c x a) 


004) 

where we have omitted s ft{ ■} since Eq. (1 1 04b is an Hermitian 
quadratic form (i.e., it is always real-valued). Therefore, the 
Gradient statistic is statistically equivalent to the GLR in 
Eq. @. 

LH statistic: In this case the general LH statistic form in 
Eq. (l65b specializes into: 


Oh = Tr [(Z wl V cA yP ai (Z wl V C}1 )Dyj 


= Tr 


z\ s-ws ,: 1 Z f 

[(PSc X a) 


Di 1 


(105) 


where D\ = Im + (Z wl V c ,iYP^ (Z wl V c ,i) in this specific 


case. Matrix D 1 can be further rewritten as: 


P»i = (l 


M 


Zt S - 1 Z e ) - 


zt STWST 3 Z, 


(a,l S c 1 a) 

Applying the Woodbury identity on -D/ 1 , we obtain: 

^r 1 ={A/ 1 + 

py (zts-'a) (zisy a )wy 


(106) 


(nt 

N (Z$Sc 1 a)tD- 1 (ZtSc 1 a)'[ 

ytl o c LI J 

a 1 Sc a 


(107) 
. Thus, 

, , - , oc»y (108) 

1 — 77 1 — 77 

with 7/ given by Eq. ( [99b . Thus LH statistic is statistically 
equivalent to the GLR for range-spread targets. 


where we exploited the definition of Dq in Eq. 
after substitution into Eq. (1 1 05b . we obtain 

2 


' M 

0 h = V + v ' 


V 


E. Standard GMANOVA 


In the present case no interference is present (t = 0, 
thus J = r). This reduces to the standard adaptive detection 
problem via the GMANOVA model considered in fTBIl . |20|| , 
ED and whose canonical form is: 


Ho: Z = N 
Hi: Z = ABC + N 


( 109 ) 
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Clearly, under the above assumptions, it holds Pa 0 = Pa 0 = 
O.Vx ,v- Therefore, the ML covariance estimate under Ho 
simplifies into Rq = K~ 1 Sq, where Sq = ZZ\ cf. Eq. (f23l> ). 

GLR: Direct specialization of Eq. ( [T 8 l > gives the explicit 
statistic: 

^ _ det[JM + (Z lvl V Cil ) 1r (Z w1 V c a)\ fllO') 

glr “ det[/ M + (Z wl V ca )t P^Z w1 V c a)} 

which coincides with the classical expressiorfj of GLR ob¬ 
tained in ca. 

Rao/Durbin statistic: Direct particularization of Eq. ( 136b 
gives: 

frao = ^{ZwqPa-l ZwoPct\ (HI) 

which provides the result obtained in ETl . 

Wald/2S-GLR statistic: Direct specialization of Eq. ( l44l > 
gives leads to: 

Lvald = T T [ZwiPA 1 Z W lPct] (H2) 


which is the same result obtained in ED. 

Gradient statistic: Direct application of Eq. (l50b provides: 


fgrad 


= 5ft Tr 


Z^PaAS^Ro 1 ' 2 : 


ZwqPc^ 


(113) 


LH statistic: In this case the LH statistic specializes into: 


t lh = Tr[(Z wl V Cil )'PA 1 (ZwiV cA )Di 1 ] . (114) 

where D\ is defined as in Sec. IIII-AI Eq. ( II 14b clearly 
coincides with the statistic obtained in |[l5l pag. 37]. 


V. Simulation results 

In Fig. Q] we report Pd vs. the SINR for all the considered 
detectors, defined as p = Tr [B^ R^^B], that is, the trace of 
the induced maximal invariant (cf. Sec. III-Ab . We underline 
that such term is also proportional to the non-centrality pa¬ 
rameter a = (0 r ,i - e rfi ) T {[i- 1 {e 0 )] g ^ 9r }- 1 (e rA - e rfi ), 
representing the synthetic parameter on which the asymptotic 
performances of all the considered test depend 0 . The curves 
have been obtained via standard Monte Carlo counting tech¬ 
niques. More specifically, the thresholds necessary to ensure 
a preassigned value of Pf a have been evaluated exploiting 
100 JPfa independent trials, while the Pd values are estimated 
over 5 • 10 3 independent trials. As to the disturbance, it is 
modeled as an exponentially-correlated Gaussian vector with 
covariance matrix (in canonical form) R = 
where er 2 > 0 is the thermal noise power, a 2 > 0 is the 
clutter power, and the th element of R, is given by 
0.95l z—J I. The clutter-to-noise ratio a 2 /<r 2 is set here to 30 dB, 
with cr 2 = 1. We point out that the specific value of the 
deterministic interference B t does not need to be specified 
at each trial considered (for both Pf a and Pd evaluation); the 
reason is that the performance of each detector depends on 
the unknown parameters solely through the induced maximal 

12 We point out that Eq. fTTob can be also re-arranged in a similar form 
as Eq. ED (i.e., a Wilks’ Lambda statistic form). Such expression, being 
equal to t g[r = det[I M + D^ 1 ^ 2 (Z w iV^,i)+ Pa 1 (ZwiV C! i) D^ 1 ^ 2 }. 
represents the alternative GLR form obtained in 03 


invariant, which is independent on B t (cf. Sec. Ill-Al l. Finally, 
all the numerical examples assume Pf a = 10 -4 . 

In order to average the performance of Pd with respect to 
B, for each independent trial we generate the signal matrix as 
B = asBg, where B g ~ CAf(O rxM ,lM,Ir ) and cub € R. 
The latter coefficient is a scaling factor used to achieve the 

desired SINR value, that is, a.b = \JP/ Tr[i?| R^\ B g \. 

For our simulation^ we assume M = 3 (i.e., an extended 
target), N = 8 , and two different scenarios of signal and 
interference lying in a vector subspace, that is (i) r = 2 
and t = 4 (sub-plots (a) and ( b )) and (**) r = 4 and t = 2 
(sub-plots (c) and ( d )). Additionally, for each of these setups, 
the cases corresponding to K = 12 and K = 19 columns 
for Z have been considered, representing two extreme case- 
studies. Indeed, the first case clearly corresponds to a so-called 
sample-starved scenario (i.e. the number of signal-free data 
required to achieve a consistent (invertible) estimate of R is 
just satisfied, that is, (K — M) = 9) while the second case 
to a setup where an adequate number of samples needed to 
obtain an accurate estimate for R is provided (i.e., in this 
case (K — M ) = 2 N = 16, with a consequent loss of 3dB 
in estimating R with the sample covariance approach, with 
respect to the known covariance case, as dictated from lf33ll ). 

The following observations can be made from inspection 
of the results. First, as K grows large, all the considered 
detectors converge to the same performance, corresponding 
to the non-adaptive case. Differently, in the sample-starved 
case (viz. the difference K — M is close to N) a significant 
difference in detection performance can be observed among 
them. First of all, the GLR has the best performance in the 
medium-high SNR range. Differently, the Rao and Gradient 
tests perform significantly better than Wald and LH tests for 
a moderate number of K (i.e., K = 12) in the case r > t 
(cf. sub-plot (c), corresponding to r = 4 and t = 2). On 
the other hand, for the same case K = 12, but r = 2 and 
t = 4, Wald and LH tests outperform Rao and Gradient tests 
when p is higher than ft! 18 dB. This is easily explained 
since both Wald (viz. 2S-GLR) and LH test both rely on 
an accurate estimate of true covariance R based on the sole 
signal-free data (cf. Secs. IIII-CllIII-FI and IIII-GI respectively). 
Differently, both Gradient and Rao tests employ a covariance 
estimate under the hypothesis Ho (that is, Rq). The latter 
covariance estimate also relies on the use of the additional 
contributions of Z corrupted by the signal 13. Although using 
them to evaluate Rq may be detrimental when the number 
of signal-free samples is adequate or the SINR is high (cf. 
sub-plots (a) and ( 6 )), when the SINR is low (i.e., the energy 
spread among the different columns is not so high) and the 
number of signal-free samples is not sufficient to guarantee 

lj Of course, due to the high number of setup parameters involved in 
the detection problem (i.e., N, K , M, r, t), we do not claim the following 
conclusions to be general for any type of setup. Nonetheless, we illustrate 
a generic setup in order to show some common trends observed among 
the detectors. A general numerical comparison is omitted due to the lack 
of space and since performance comparison in some specific setups (such 
as those considered in Sec. |IVt can be found in the related literature. 
Nonetheless, the supplementary material attached to this paper contains some 
additional numerical results aimed at confirming the statistical equivalence 
results obtained for the considered special scenarios. 
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(a) (b) 



Figure 1. Pd vs. p for all the considered detectors; common parameters: M = 3 and N = 8. Case (a) (top-left) r = 2, t = 4 and K = 12; Case ( b ) 
(top-right) r = 2, t = 4 and K = 19; Case (c) (bottom-left) r = 4, t = 2 and iC = 12; Case (d) (bottom-right) r = 4, t = 2 and iC = 19. 


a reliable estimate of R, the degradation of using signal- 
corrupted terms is overcome by the (beneficial) availability 
of additional samples for covariance estimation. 

VI. Conclusions 

In this second part of this work, we have derived several 
detectors for adaptive detection in a GMANOVA signal model 
with structured interference (viz. I-GMANOVA). We derived 
the GLR, Rao, Wald, 2S-GLR, Durbin, Gradient, and LH 
statistics. All the aforementioned statistics have been shown to 
be CFAR with respect to the nuisance parameters, by proving 
that all can be written in terms of the MIS (obtained in the 
first part of this work). For the considered general model, we 
also established statistical equivalence between: ( i ) Wald and 
2S-GLR statistics and ( ii ) Durbin and Rao statistics. 

Furthermore, the following statistical-equivalence results 
have been proved in the following special setups: 

• For point-like targets (with possible point-like interfer¬ 
ence), we have shown that Gradient and LH tests are 
statistically equivalent to Kelly’s GLRT; 

• For multidimensional signals, we have shown that: (a) 
Rao test is statistically equivalent to Gradient test and 
( b ) Wald test (2S-GLRT) is statistically equivalent to LH 
test; 


• For range-spread targets and rank-one subspace (r = 1), 
we have shown that Gradient and LH tests are statistically 
equivalent to the GLRT. 

Finally, simulation results were provided to compare the 
performance of the aforementioned detectors. 

VII. Supplementary Material Organization 

The following additional sections contain supplemental ma¬ 
terial for part II of this work. More specifically. Sec. IVIIII 
contains the proof of Lem. [j] in the paper, while Secs. IIXI 
IXl and IXll provide the derivation of Rao, Wald, and Gradient 
(Terrell) tests, respectively. Furthermore, Sec. IXIIl provides the 
statistical equivalence between Rao and Durbin tests (Thm.|2]i. 
Additionally, Sec. IXIIll nrovides a series of useful equalities for 
showing the CFARness of all the considered detectors. Finally, 
Sec. IXIVl provides some numerical results aimedat confirming 
the special equivalence results obtained in the manuscript. 

VIII. Proof of Lemma 1 

We only provide the proof for Eq. (l24l >. as the equality for 
R^Et in Eq. ( l25l > can be obtained following similar steps. 
We first rewrite Eq. d22l> as: 

Rr = K~ l [S c + Sl /2 [PX 1 Z w 1 )Pc\{Pa 1 Z w PS\ t 2 } 

(115) 
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Taking the inverse and exploiting Woodbury identity 
gives: 

Ry 1 = K [S- 1 - S-^iP^Z^Vd 
x{I M + {Z W1 V C ^ P^iZ^V^)}- 1 

xV^Zl.P^S- 1 ' 2 ] (116) 

It is apparent that the second term in Eq. (II 1 6b is null when 
post-multiplied by A , since (P^Sc 12 A) = P A A\ = 
Onxj, thus leading to the claimed result. 


IX. Derivation of Rao statistic 
In this appendix we report the derivation for Rao statistic in 
Eq. (SB . Before proceeding, we define the auxiliary notation 

b S M = [bn bf Ji ] T and b sJ = [bj b 
shown that: 

din MZ: B S ,R) 


'll] 


. First, it can be 


db StR 

din 

db a _T 


25ft {g A } 
2 5ft{gs} 

2 3{9a} 
2 3{gs} 


where: 


£ R JMx1 ; 

(117) 

£ R JMx1 ; 

(118) 

irM X1. 

(119) 

itMxl 

(120) 


g A = vec (El R 1 Z d C7+) £ C r 

g B = vec (El R 1 Z d C*) £ C‘ 

In Eqs. (E3 and (fHOt . we have adopted the simplified 
notation Z d = (Z — AB S C). The results in Eqs. (1117b and 
(II 18b are obtained by exploiting the following steps: 


1) Evaluate the complex derivatives 91n 
and o^mz^.r) (as wdl as am Mz-.b.,r) 


and 


dinf^z-^B^R) ) stan d ar d complex differentiation 

rules 041 : 

2) Exploit that for any f(x) : C pxl —> R, it holds 
SM_= and Kg) = 23{«g>( (see 


e-g- 


351); 


3) Obtain ,s u 


and 


a In h(Z-,B 3 .R) 


[ s,R db B j 

tion of the gradients obtained at step 2). 


as composi- 


By identical steps, it can be also proved that: 


d\nfy(Z-B s ,R) din fy(Z-O) 


dO r 


dO r 


25ft {g A } 
2 


( 121 ) 


which gives the explicit expression for the gradient of the log- 
likelihood required for evaluation of Rao statistic (cf. Eq. (l29l)). 
On the other hand, the block (0 r ,0 r ) of the inverse of the FIM 
is evaluated as follows. First, we notice that 051 : 

= l/o 1 ( 0 )]g r ,9 r d22) 

where [ I a (6 )] is here used to denote the block of the 
FIM comprising only the contributions related to (0 r ,0 St a). 
The aforementioned property follows from the cross-terms 
( 9 r ,6 s .b ) and (0 s , a , Os,b) being null in the (overall) FIM 
1(0). Additionally, the following equality holds (recalling 


(flTTb and fel l: 


that 1 01 OU = 


[b T n 




bln 


bf j 1 and exploiting Eqs. 


8 In MZ;B S ,R) 

d 

' 0 r ' 

0s,a 



= P 


din f 1 (Z:B s ,R) 

d In fi(Z-,B s ,R) 
db s t 


(123) 


where P £ ]|2 JMx 2JM j s a su jt a bl e permutation matriJT 
defined as 


P = 


IrM 

OrMxrM 
0 tMxrM 
0 tMxrM 


QrMxtM 

OrMxlM 

ItM 

0 tMxtM 


OrM xrM 

IrM 

OiM XrM 
OiM XrM 


0 rMxtM 
OrMxiM 
0 tMxtM 
ItM 


(124) 


Before proceeding, we define the matrix fi = (A^ R 1 A) and 
the partitioning: 


fill 

fil2 


~E\R- [ E t E\R- l Er 

fi21 

fi22 


ElR-'Et ElR- 1 E r _ 


(125) 


where, fi,y, ( i,j) £ {1,2} x {1,2}, is a sub-matrix whose 
dimensions can be obtained replacing 1 and 2 with t and r, 
respectively. Then, the sub-FIM I a (0) is obtained starting 
from Eq. (11231) as I a (0) = (Pf P T ), where \1/ has the 
following special structure: 


* = 


2 5 R{K} 

2S{1T} 


-2 %{K} 

2 5 ft{K} J ’ 


(126) 


and the matrix K £ C JMxJM j s defined as: 




(cc t) T ® fi 22 
(cct) r ®o 12 


(cci) T ® n 2 i 
(cci) T ® 


(127) 


Finally, the inverse I~ x (0) is obtained as I~ x (0) = 

(P\&“ 1 P T ) (as P is orthogonal). In the latter case, the 
inverse matrix \£'~ 1 has the same structure as fir (cf. Eq. ( 1 126b ) 
except for K and the factor 2 replaced by K" 1 and }, 
respectiveh{3 By exploiting the following block structure of 
K l 


K 1 


K n 

K 

K 21 

K 


(128) 


with K 11 G IT 12 G C r MxtM 21 ^ ^ tMxrM 

and K 22 £ C tMxtM , respectively, and the structure of P, it 
can be shown that: 




|s{K n } ^{K 11 } 


(129) 


Similarly, it is not difficult to show that K 11 is given in closed- 
form as: 


k 11 = {(cci) T ® fi 22 - {ccftf ® n 2 i 

X [(CC^) T ® fin] ~ 1 (CC ] ) T <3 fii 2 } (130) 
= (cc t )- T ®r 22 (131) 

where Igy is a sub-matrix obtained from T = fi -1 exploiting 
identical partitioning (in terms of size) as done in Eq. (1125b for 
fi. The compact expression in Eq. (1131b is obtained from the 
use of mixed product and associative properties of Kronecker 


14 Recall that every permutation matrix is a special orthogonal matrix, that 
is, P 1 = P T . 

15 Such result is obtained by exploiting the equality V F 1V I' = I and the 
real-imaginary parts decompositions of V F = 'I' /,■ + j V I' / and V F 1 — '1' /,> ■ 
j T 7 /, from which it follows the set of equations (i) ( y F i; y F /,■ — d/ / v l/ /) = I 
and (ii) = 0. 
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operato0 Then, combining Eqs. (1121b . ( 1 1 29b and ( 11311 ) leads 
to: 


j ? In fi(Z;0) 


w], 


a In /i(Z;0) 


dQT r ^ne r ,e r QQ r 

vec(El R~ x Z d C/ [(CC t )- T ®r 2 


ex 


xvec 


(£?t JJ-iZjCt) 


= Tr 


Z\ R 1 E r T-n E'lR- 1 Z,iP, 


d-^Ci 


(132) 

(133) 


where we have exploited the well-known equivalence between 
a real-valued Hermitian quadratic form and its real symmetric 
quadratic counterpart in Eq. (11321 ) and some standard proper¬ 
ties of vec(-) operator^ in obtaining Eq. (11331 ). Finally, the 
substitution 0 = 6q provides: 


Therefore, in view of Eq. (11391 ). it holds: 


p i o R- l/ 2 E r T? 22 ElR- il 2 Pi o 
(R/ /2 E t )(T ° 21 ) f + (. R~ 1/2 E r ) (f° 2 )tl (f° 2 )' 


(140) 


r°AEj R- 1/2 ) + T° 22 (ElR^ 1/z ) 
= Rq 1 / 2 A(A^Rq 1 A)- 1 A^R / /2 


i?t R7 1/2 ' 


(141) 


-(Ro 1/2 E t )(ElR^E t )-\ElR ~ 1/2 ) (142) 

= (Pa 1 ~Pa 0 ) d43) 

where we have exploited the equality {E\r/ E t y A = T^ — 
Ti2(T22) _1 T 21 (which can be deduced from Eq. (11251 ) after 
substitution R = Rq and from T° definition). Finally, in 
Eq. (11431 ) we have further defined A\ = (R a l ^ 2 A). 


Tr 


z\ fi R/E r T° 2 


P'1 R q 1 Z d ,0 Pci 


(134) 


where Rq and Z d ,o are defined in Eqs. ( l23l > and (l33l >. 
respectively. Similarly, TT denotes a sub-matrix obtained from 
T°j = (A'R {) 1 A) -1 by exploiting identical partitioning (in 
terms of size) as done in Eq. (11251 ). This provides the explicit 
expression for the Rao statistic. 


Proof of Eq. ( P?l ) 

The mentioned result is proved as: 

R~ 1/2 Z dfi P c t 

= Rq 1 ' 2 (Z - Si / 2 P Ao 5-V 2 ZP ct ) P c t (135) 
= P~ 1/2 (l N - Sj/ 2 P Ao S - 1 / 2 ) Z P ct (136) 

where we have exploited Eq. (1331 and Pci Pci = Pet, 
respectively. The above expression can be further rewritten 
by exploiting Eq. (l25l > of Lem. |T| as: 

R- ll2 Z dfi Pa 

= (Rf 1/2 - R/^EtiEjs/E^E} S/)ZP C t (137) 

= ( R / 12 - P Ao R/ /2 )Z P C i = P/R/ /2 ZP C i 

(138) 

where we have denoted To = R Q 1 _ E t , which finally 
provides Eq. (l34l > (after having defined Z wo = R 0 1 2 Z). 


Proof of Eq. ( 071 ) 

First, it can be shown that: 


T° 2 El R/ /2 Pi = (f° 21 E\ + f ° 2 El) R/ /2 , 


039) 


X. Derivation of Wald statistic 
In this section we report the derivation for Wald statistic 
in Eq. (l42l >. In order to prove the aforementioned result, 
we build upon the explicit expression of [l~ 1 (0)] e e ob¬ 
tained in Eq. (11291 ). Such result allows to readily evaluate 
{[ l~ 1 (0\)]e r .e r }~ 1 in Eq. SB by ( i ) substitution R = Ri 
and ( ii ) matrix inversion 1 1 as: 

WW 1 )- 1 } -2SKX1 1 )- 1 }' 
^{(iTj 11 )- 1 } 25R{(K 1 11 )- 1 } J ’ 

(144) 


{[T- 1 ^ i)] 0 rt e r }- 1 = 


where (K\ x ) 1 = (CCT) T ® (T 22 ) 1 , and TT is a sub¬ 
matrix obtained from T 1 = (A^ R/ A/ 1 exploiting identical 
partitioning (in terms of size) as done in Eq. (11251 ) for D. 

We have now to evaluate Q r o and 0, i, respectively. 
First, we recall that 6 


r, 0 


= o 


n T 


2rM, 


while d r 


5i{vec(P)} T T{vec (B)} T , with B representing the ML 

estimate of the complex-valued signal matrix under Hi. This 
estimate can be shown to be equal to: 

B = K r 22 El S/ /2 P/S/ /2 ZC f (CCl)- 1 . (145) 

The above result is obtained starting from Eq. (fl4l) and 

n T 


observing that B s = 


BTi 


B 1 


. Therefore, collecting the 


above results, Wald statistic is obtained as 

(Or, 1 - Or,o) T (Or, 1 - Or , 0 ) 

= 2 vec(P)l (iTj 11 ) -1 vec(B) (146) 


oc vec 


(?* 2 El S// 2 P/S/' 2 Z Ct(CCt) 


-1 


xK {(CCt) T ®(fi 2 )- 1 } 


xvec 


= Tr 


(147) 


(fi 2 El S/I 2 p/S/I 2 zc\c&y x 

Z/(KP/s/> 2 E r f\ 2 ElS// 2 P/)Z wl Pci^ 

(148) 


16 The mixed product property states that (V 10 V 2 XV 30 V 4 ) = (Vi Vs)0 
(V 2 V 4 ), where Vi are generic matrices of compatible sizes. Differently, the 
associative property states that V± 0 (V 2 + V 3 ) = Vi 0 V 2 + Vi 0 V 3 . 

17 More specifically, we have exploited vecfViV^Vs) = (V^ 0 

Vi) vec(V 3 ) and vec(Vi)^ vec(V^) = TrfV^ V 2 ], with Vi being generic 
matrices. 


where we have again exploited (as in the derivation of Rao 
statistic) the well-known equivalence between a real-valued 

18 We again use t he pro perty that the inverse of a block-symmetric matrix 
in the form of Eq. n~29i gives rise to a similar structure for its inverse, as 
exploited for the derivation of Rao statistic. 
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Hermitian quadratic form and its real symmetric quadratic 
counterpart in Eq. (1 1 46b and some standard properties of vec(-) 
operator in obtaining Eq. ( 1148b . This provides the closed form 
expression for Wald statistic. 


Proof of Eq. ( I?i1 ) 

We first notice that the following equality holds: 

f\ 2 El S " 1 / 2 Pi 0 = (f h E\ + f\ 2 El) Sf 1 ' 2 . (149) 

The above result is an almost evident consequence of the 
application of matrix inversion formula for a 2 x 2 block 
matrix. Therefore, in view of the above equality, we observe 
that: 


K P^Sf 1 / 2 E r f} 22 ElSf^Pl 
= K [(S- 1 / 2 ^)^!) 1 + (S-^ 2 E r ) (f£ 2 )t' 

x [rl^ElSf 1 ' 2 ) + T 1 22 (ElSf 1 / 2 ) 
= ( Pa i - Pa 0 ) = Pa 
which thus provides Eq. ( l43l >. 


(r^r 1 


(150) 

(151) 


XI. Derivation of Gradient statistic 


which readily follows from application of Eq. ( l34b . Then, we 
rewrite the matrix R Q P^R 0 1 as: 


D — 1/2 rp_L d 1/2 

= Sf 1/2 [sl^Ro 1 - Sy 2 Rq 1/2 P Ao Rq 1/2 
= S-W [l N - Sy 2 {R^E t ) 
x {ElR^E^ElSf 1 ' 2 ] S l J 2 Rf l 

— Sf 1 / 2 p Ao s l J 2 Rf 1 


(156) 


(157) 

(158) 


where in Eq. ( 1158b we have exploited Lem.Q] Eq. ( l25b . Finally, 
from straightforward combination of the results in Eqs. (1155b 
and (1158b . the final result follows. 


XII. Equivalence between Rao and Durbin 
statistics (Theorem [2]) 

In this section we prove statistical equivalence between 
Rao and Durbin statistics by explicitly deriving the closed 
form expression of Durbin statistic (implicitly expressed in 
Eq. (l55b ). To this end, analogously as for the Rao and Wald 
statistics, it can be shown that: 


The derivation of Gradient statistic is readily obtained from 
intermediate results obtained in derivation of Rao and Wald 
statistics, in Secs. [IX] and [X] respectively. Indeed, exploiting 
Eqs. (1121b and (1 1 45b provides: 


din h(Z-O) 


dOj 

oc iGR 1 TV 


(flr,i - 0r,o) = 2K{vec(B) t g^} (152) 


0=0o 


(fi 2 El S " 1/2 Pi 0 Sf^ 2 ZC\CCy~ l ) ] 

xElRf' Z dfi C 1 ] } (153) 

= 3? {Tr [zl t KPi a Sy/ 2 E r f\ 2 ElKf 1 Z d , 0 P ct ] } 


(154) 


In Eq. ( 1152b we have exploited the equivalence between the 
real part of an inner product in the complex domain and its 
real-valued equivalent counterpart*^. along with the defini¬ 
tion g° A = vec(El Rf 1 Zd t o C^). Furthermore, in obtaining 
Eq. (Tf53l we have exploited standard properties of vec(-) 
operator. Finally, we recall that Zdfi is given in Eq. ( l33b . This 
concludes the derivation of Gradient statistic. 


Proof of Eq. d?9l > 

We start by observing that 

Rf 1 Z d fiP c , = Rf 1/2 Py o Rf 1/2 Z P ct , (155) 

19 More specifically, given two complex vectors with real/imaginary parts 
decomposition v\ = v\ # + jv i j an 112 = v 2 R 3 V 2 /> it holds 

5R{uji>2} = 5R{u|ui} = vf e V2,e, where v itE — [ v T,r v T,i\ ■ 



where: 


b K{t 0 } 

.1 Z{T 0 } 
2 5R{T 0 } 
2 4{T 0 } 


-h*{T 0 y 

i5?{To} 
-2 3{T 0 }' 
2 5R{T 0 } 


059) 

060) 


To4(CCl)- T ®f° 2 ; (161) 

To = {CC ] ) T ® {El R 0 1 E r ). (162) 


The result in Eq. ( 1159b is obtained starting from the explicit 
expression of [/ -1 (0)] in Eq. ( 1129b and plugging back 

R = R<). Differently, the estimate 0 r ,oi can be obtained 
following the steps described next. Without loss of generality, 
we consider maximization of ln(-) of the objective in Eq. (f56t 
and evaluate the following estimate: 


Bn = 


(163) 


arg min Tr 

B 


Z — A 


B t .o 

B 


C) Rr\ 


Z A 


Bt.o 

B 


Once we have obtained Bq, 0, oi is evaluated through the 

r ’ _ A iT 

simple operation 0 r , O i = 3ff{vec(JB 0 )} T 3{vec(-B 0 )} T . 

Thus, it is not difficult to show that the solution to the 
optimization problem in Eq. (1163b is given in closed-form as: 


B, 


= (ElRo 1 


El R^ 1 Z, 


< 2,0 


C*(CC* 


(164) 
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Substituting Eqs. d 1 59b , ( 11601 ) and d 1 64b into Eq. d55K pro¬ 
vide© 


oc vec 


it 


xvec 


2 vec(-B 0 ) t (T 0 T 0 T 0 ) vec(^ 0 ) 

(#^#o 1 Ery 1 El % 1 Z^oC'iCC')- 1 
{(CCf <® [(^tHo ” 1 E r ) f ° 2 (^eIRq 1 JS r )] } 
(eIr^ 1 E r )El Ho 1 Z d< o C\CC 




= Tr [Z\ o Hq 1 H r f ° 2 El Rf 1 Z d , 0 P« 


c H 


(165) 

(166) 


where we have exploited standard properties of vec(-) in 
Eq. (fT66t . 

Clearly, the last result coincides with the Rao Test for the 
I-GMANOVA model, as apparent from comparison with Rao 
statistic reported in Eq. ( |3H in the manuscript. 


XIII. Proofs of useful equalities 
Proof of Eqs. ( 1261 ) and ( 1271 ) 

Hereinafter we provide the proof of Eqs. ( l26t and < [27l >. 
which are fundamental for proving CFARness of GLR, Wald 
and LH statistics. We start by observing that: 

(Z wl V Ctl )^ P^ 0 (Z W1 V C<1 ) 

= Z\ {s- 1 - {S-^E t ){E\s- c x E t )~ x E\ S" 1 } Z c (167) 
Analogously, we write: 

(Z wl V Cll )i Pjt^ZnVc, t) 

= Z\ {5- 1 - (S^A^S^A)- 1 ^ S" 1 } Z c (168) 


Before proceeding further, we define the following partitioning 
for matrix S ’" 1 as: 


'll 

S 12 

s 13 

'21 

S 22 

s 23 

'31 

s 32 

s 33 


(169) 


Furthermore, S'f (i,j) € {1, 2, 3} x {1, 2,3}, is a sub-matrix 
whose dimensions can be obtained replacing 1, 2 and 3 with 
t, r and (N — J), respectively. First, exploiting E t structure 
leads to: 


(EjSf 1 E t ) = S 11 , 


(S-'Et) 


'S 11 ' 
S 21 
S 31 


(170) 


Accordingly, the matrix within curly brackets in Eq. (1 1 67b can 
be rewritten as: 


{ 


7-1 


(S^E^El 

0 txt 

P(N-t)xt 


s-'E t )-'E!s;f 

o tx(N-t) 

0-1 ) 

°2 J 


(171) 


(172) 


where we have denoted: 

^ A_ S22 S 23 

2 “ [S 3 2 5 33 . ' 

Secondly, after substitution in Eq. < 1 1 67b . we obtain: 

(Z W iV c ,i)i P^ 0 (Z wl V cA ) = Z \ 3 Sf 1 Z 23 ; (173) 

where Z 23 = [Zf Z : 'f ] T . Finally, by exploiting the block 
inverse expression of S 2 in Eq. (1173b . it follows that 


(Z wl V cA ) f Pi 0 (Z wl V cA ) = zl 3 Sfl Z 2 .3 


A Sf 3 l Z 3 


, 

(174) 

which proves Eq. ( l26l >. Similarly, exploiting A structure, it can 
be shown that: 


A*S7 1 A = 


S 11 

s 21 


s 12 ' 

S 22 


S7 1 A = 


' s 11 

S 12 ' 

s 21 

s 22 

_s 31 

s 32 _ 

within 

the 


(175) 


brackets in Eq. ( 1168b as: : 

{S " 1 - (S^A^S^A)- 1 ^ S 7 1 } 

0 JxJ 0 Jx(N-J) 


( N-J)xJ 


c-i 
^33 


Finally, gathering the above results leads to: 

(z wl v Ctl y Pi^z^v^) - ^ 

which proves Eq. (l27l ). 


A S£ z 3 


(176) 


(177) 


Proof of Eqs. ( li<SI ) and ( U9I ) 

Hereinafter we provide a proof of Eqs. (l38l > and ( l39l ). 
which are fundamental for proving CFARness of Rao (Durbin) 
statistic. Firstly, it can be shown that: 


D—1/2 nl p—1/2 

■^o Ft„ u o 


S—1/2 pi S 1/2 _ 

#0 “a 1 #0 — 


vtxt 

.0 (N-t)xt 

0 JxJ 

P(N-J)x J 


0 


tx(N-t) 

# 


■0,2 


0 


Jx(N-J) 


H, 


■0,33 J 


where we have defined the following partitioning: 


(178) 

(179) 


Ho. n 

Ho, 12 

Ho, 13 

Hq. 21 

Ho, 22 

Ho, 23 

#0.31 

Ho, 32 

Ho, 33 


Ho. 22 

Ho, 23 

Ho. 32 

#0,33 


(180) 


where Ro,ij, (i,j) 6 ( 1 , 2 ,3} x {1,2,3}, is a sub-matrix 
whose dimensions can be obtained replacing 1, 2 and 3 with t, 
r and (N — J ), respectively. Then, exploiting the above results 
we obtain: 


(Z W 0 V Ctl )i Pj[ o (Z wo V C!l ) = Z 2 t 3 H 0 -^Z 23 (181) 
(Z W 0 V r c , 1 ) t (Z W 0 V cA ) = Z\ H^gg Z 3 (182) 

Furthermore, it is not difficult to show, starting from Eq. ( |23| >. 
that 


20 We have exploited the fact that To To To is Hermitian and that the 
product of block-symmetric real counterparts of Hermitian matrices can 
be expressed as an equivalent block-symmetric real counterpart with the 
component matrix being given by the product of the aforementioned matrices. 
Finally, we have used the equivalence between an Hermitian quadratic form 
and its real block-symmetric counterpart. 


Ho ,2 = K - 1 \S 2 + Z 23 Zl 3 \, (183) 

# 0,33 = K - 1 [S 33 + Z 3 zl\. (184) 

Finally, after substitution into Eqs. ( 1181b and (1 1 82b and appli¬ 
cation of Woodbury identity, the claimed result is obtained. 
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Proof of Eqs. ( 1521 ) and ( I5JI ) 

In what follows, we provide a proof of Eqs. (l52l > and (l53l >. 
which are exploited in Sec. I1I1-DI for proving CFARness of 
Gradient statistic. We first observe that the following equalities 
hold: 

(Z W1 V C ^ Pi 0 S p2 R~ 1/2 (Z W0 V Ctl ) 

= {Z W0 V cA y Pi 0 {Z W0 V c>1 ) (185) 

(Z wl V cA y Pf S 1 / 2 Rf 1/2 (Z W0 V cA ) 

= K {{Z wl V c ^Pf (Z W1 V C>1 ) 

-(ZnVcrfP^ Z wl V cA (. I M + (Z wl V cA ) f Pi 0 Z wl V cA ) 
x(Z wl V Ctl ^Pi o Z wl V cA ] (186) 

where we have used Eq. (l25l > of Lem. [T| in deriving the right- 
hand side of Eq. (1 1 85b . Differently, Eq. < 1 1 86b is obtained ex¬ 
ploiting the following steps: (i) use of R<, explicit expression 
given by Eq. ( |23| >. (it) application of Woodbury Identity to 
Rf 1 and (in) simplification through the use of the equality 
Pj^ t Pa 0 = p a, ■ Finally, by exploiting Eq. < 13 8 b into Eq. (11851 ) 
and Eqs. (l26l> and (127l) into Eq. (1186b . we demonstrate the 
considered equalities. 

XIV. Simulation results showing specific 

COINCIDENCE RESULTS 

In this section we confirm, through numerical results, the 
statistical equivalence results obtained among the considered 
detectors for specific adaptive detection scenarios. We remark 
that the simulation parameters (i.e., the the structure of the 
covariance R and the generation process for unknown signal 
matrix B ) are the same as those used in the manuscript (as 
well as the Monte Carlo setup) and thus are not reported for 
the sake of brevity. 

First, in Fig. U we show Pd vs. p (given Pf a = 10~ 4 ) 
for a setup with point-like signal and interference (M = 1) 
where the signal belongs to a two-dimensional vector suspace 
(r = 2) while the interference to a four-dimensional vector 
subspace (t = 4). We assume that each column of Z is a 
vector of N = 8 elements and K = 13 samples are assumed. 
It is apparent the statistical equivalence among GLR, Gradient 
and LH tests. 

Similarly, in Fig. [3] we show Pd vs. p (given Pf a = 10~ 4 ) 
for a setup with multidimensional signals (N = r = 8) where 
M = 8 and K = 24. It is apparent the statistical equivalence 
between Rao and Gradient tests and between Wald and LH 
tests. 

Finally, in Fig. [4] we report Pd vs. p (given Pf a = ICC 4 ) 
for a range-spread target (M = 8) with a rank-one signal 
subspace (r = 1) and no interference (t = 0), where K = 24 
and N = 8. It is apparent the statistical equivalence among 
GLR, Gradient and LH tests. 
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